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1. Introduction 



A Probabilistic Cellular Automaton (PCA) [MS1,LMS,GD] , or interacting par- 
ticle system [Lig], can be regarded as an infinite collection of cells or particles, 
where each cell or particle can take on a finite number of states. The discrete time 
evolution of such systems is determined by a statistical law according to which, in 
any given configuration at time t, all cells are simultaneously and independently 
updated to the configuration at time t + 1. One of the basic questions concerning 
such systems is ergodicity, i.e., the uniqueness of the stationary probability measure 
on configurations. In this article we will introduce a quantum analogue of this struc- 
ture, and we will also prove an analogue of well-known criteria by which ergodicity 
of a PCA can be decided in terms of the local transition probabilities [Lig,LMS, 
MSI]. 

Our definition of quantum cellular automata (QCA) is an abstraction of a 
structure that arose in the project of generalizing the construction of "finitely corre- 
lated states" on quantum spin chains [FNW1,FNW2] to two-dimensional systems 
[Ric]. In this context the ergodicity of the QCA is equivalent to the state on the 
two-dimensional system to be uniquely determined by local data, independently of 
boundary conditions. The term "quantum cellular automaton" has been used pre- 
viously by some other authors [GZ,FGSS,LTPH] . In the cases we are aware of, 
however, it is used for a structure on the Hilbert space level, and not on the level 
of observables. Thus in [GZ] the classical states at each site are simply replaced 
by the values of the wave function at that site, and the dynamics is just a discrete 
Schrodinger equation with non-selfadjoint Hamiltonian, made non-linear by keep- 
ing the normalization fixed. There is some interest in quantum cellular automata 
also from the point of view of nanometer scale computers, for which quantum ef- 
fects are expected to be relevant [Mai,Bia,LTPH]. The evolution of the automata 
considered in this paper is in general non- unitary, i.e., pure states may evolve into 
mixed states. This might be an interesting addition to the structure of "quantum 
computers", as studied by a number of authors recently (see [DiV], and references 
cited there). 

The main technical contribution of this article is the introduction of a general 
class of "oscillation norms" on quantum lattice systems. We believe this to be 
a useful tool of independent interest, and therefore included proofs of the basic 
general properties of such norms. They generalize a classical notion, which also 
was a principal tool for the proof of the ergodicity criteria in [LMS,MS1,MS2]. 
A special case of such a norm in the quantum case was already used extensively 
by Matsui [Mal,Ma2,Ma3,Ma4]. Among other things, he used it to establish 
ergodicity criteria for the continuous time analogue of QCAs. Another special case 
of oscillation norms, used for the same purpose, is to be found in a recent preprint 
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by Majewski and Zegarlinski [MZ]. 

An important consideration in the study of transition operators on composite 
quantum systems is that the norm of an operator acting on observables may increase, 
if we consider the given system as a subsystem of a larger one. It is therefore es- 
sential to consider versions the basic operator properties, which are "stabilized with 
respect to system enlargement" . The stabilized versions of positivity and bounded- 
ness (called "complete" positivity and boundedness) are well-known, and we give a 
brief summary of these with references in an Appendix. For the notion of "bound- 
edness in oscillation norm" the stabilized version is described in Section 4. Much 
to our surprise, it turned out that in this case the necessity of stabilization is not 
characteristic of the quantum case. Even in the classical case, as soon as one has 
more than two states at each site, the oscillation norm bound of a transition oper- 
ator may increase with the size of the environment, as we will show by an explicit 
example. 

The paper is organized as follows. In Section 2 we establish our notations 
for quantum lattice systems, and develop the definition of QCAs from the classical 
notion of PC As. Section 3.1 introduces oscillation norms in general C*-algebras, and 
shows how contractivity of a transition operator in such a norm entails ergodicity. 
Section 3.2 is devoted to the construction of a canonical oscillation norm for a 
system composed of many parts with given oscillation norms, such as the quasi- 
local algebra. The key result, reducing estimates on the infinite systems to an 
estimate of a local quantity, is shown in Section 4.1: if the local transition operators 
contract with a certain rate, as measured by the "completely bounded oscillation 
norm", then the same holds for an infinite tensor product of such operators. We 
briefly indicate in Section 4.2 how other plausible approaches to such estimates fail 
on this account. The basic ergodicity criterion for QCAs in Section 5.1 is a direct 
corollary of this estimate. In 5.2 we show in what sense the classical PCAs are 
covered by this criterion. Finally, we show, in Section 5.3, that the same estimate 
on local transition operators, which implies ergodicity by our general criterion, also 
entails exponential decay of correlation functions in the unique invariant state. A 
quick review of the notions of complete positivity and complete boundedness, with 
pointers to the literature, is given in the Appendix. 
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2. Definition of Quantum Cellular Automata 




In order to describe the notion of quantum cellular automata (QCA), it is best to 
begin by restating the classical structure in an algebraic language more suitable for 
generalization to the quantum case. The underlying lattice, or set of cells, will be 
denoted by C. To each cell x G £ we associate an observable algebra A x , which in the 
case of a classical system is simply the algebra A x = C (fi x ) of continuous complex 
valued functions on the set O x of configurations of each cell. Finite subsystems, 
associated with finite subsets Ac£ are described by the tensor product 



(2.1) 



This formula is also used for infinite subsystems, and in particular, for the observable 
algebra A c of the whole system. The product on the right hand side is then the 
infinite Cartesian product of compact topological spaces. Even if each cell has 
only finitely many configurations, and hence "continuity" for functions on each fi x 
is a vacuous condition, continuity of the observables / G A c is the non-trivial 
requirement that / may be uniformly approximated by observables depending only 
on finitely many cells. 

The PC A dynamics is determined by the transition functions 

p x (u,a) = Eju^t+i G a u t = u;} , 

for x G £, lo G O, and measurable a C O x . This is a probability measure on Q x 
in its second argument, depending, in principle on the entire previous configuration 
w£fl. In algebraic language this becomes an operator 

P ■ A x — > A c 

(2.2) 



(P x f)(uj) = J p x {u),du x ) f(u x ) . 



The assumption that P x f is a continuous function is called the Feller property of 
the PCA [Ma4,Lig]. It means intuitively that the updating of one cell does not 
depend too sensitively on infinitely many other cells. It is automatically satisfied 
for finite range interactions, i.e., when P X (A X ) C A A ^ C A c , for some finite set 
A(x). 

Independent updating of different cells, the basic property of PCAs, means 
in more formal language that the distribution of uJt+i, given ui t , is the product 
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measure formed out of the measures p x (u> t , •). Equivalently, the transition operator 
P : A c — > A c is defined as 




where f x G A x , A C C is an arbitrary finite set, and the product on the right 
hand side is the pointwise product of functions in the algebra A c = C (O). Since the 
tensor products on the left hand side of (2.3) span the C*-algebra C (fi), this equation 
determines P uniquely. This concludes our brief description of PCA dynamics. 

Some of the above is easily translated into the quantum setting: the main 
change is that now all observable algebras may be non-commutative (rather than 
commutative) C*-algebras with identity. For the observables algebra of a single cell 
one typically chooses the algebra M. n of n x n-matrices with n < oo. The simplest 
example is a Heisenberg spin- 1/2 system, for which A x = M.2 for every x G C. 
The observable algebra of a composite quantum system is defined as the closure 
of the algebraic tensor product in a suitable C*-norm. In contrast to the classical 
case there may be several such norms, in which case we always take the "minimal" 
C*-norm [Tak] . For the product of finite dimensional matrix algebras, in which we 
will mostly be interested, all C*-tensor norms coincide anyhow. We continue to use 
the notation A A for the observable algebras (2.1) of finite regions. For A t C A 2 , 
there is a natural inclusion A Al C A^ 2 , by tensoring each element of A Al with the 
identity in A Ai \ Al . The infinite tensor product defining the observable algebra of an 
infinite (sub-) system is always defined as the C*-inductive limit [BR] of the finite 
tensor products with respect to these inclusions. As in the classical case, this simply 
means that all observables can be approximated in norm by finitely localized ones. 
A transition operator such as P x can be characterized as a positive operator (i.e., 
an operator taking positive elements into positive elements), mapping the identity 
into the identity. Moreover, we will also assume that these properties persist if we 
consider A x and A c as subsystems of a larger system, i.e., we require P x to be 
completely positive (see the Appendix for definitions). A probability measure is 
replaced in the algebraic framework by the expectation value functional, or "state" 
it induces. Thus states are linear functionals on the observable algebra, which take 
positive values on positive elements, and the value 1 on the identity. 

The key problem for the quantum generalization of PCAs is the positivity of 
the right hand side of equation (2.3): in the non-commutative context a product of 
positive elements is practically never positive. In fact, restricting to the case of just 
two factors (A = {x, y} with f x ,f y > 0), we find that a necessary and sufficient 
condition for the positivity of P is that the ranges of the operators P x commute in 
A c (see Proposition IV.4.23(ii) in [Tak]). 

We therefore have to choose our definition in such a way that the commuting 
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range condition holds automatically. The following is one way of doing this. It is 
perhaps not the most general possibility, but it covers the cellular automata which 
came up naturally in our construction of states on two-dimensional spin systems. 
The idea is to subdivide each cell into subcells, such that the images of different 
P x are contained in different tensor factors with respect to the refined tensor de- 
composition of A c . For notational convenience we state the definition only in the 
case that the observable algebras and transition operators of different cells are all 
isomorphic. 

1 Definition. A quantum cellular automaton (QCA) is given by the following 
objects: 

(1) a countable lattice C, 

(2) a set S of "subcell types", and a C*-algebra B s with unit, for each type s G S, 

(3) for each s G S, an injective "propagation map" d(s; •):£—>£, and 

(4) a completely positive unit preserving operator Pi : <S) seS B s — > ® s£ s^ s - 
The following objects are defined in terms of the above: 

(5) at each site x G C the algebra A x = <S> seS B x ' s , where B x ' s is an isomorphic 
copy of B s , 

(6) the quasi-local algebra A c = < S> xe cA x , 

(7) the transition operators 

P x :A x = (g)& c ' a — > (g)B d(s ' x ^ s CA C , 

ses ses 

where P x is defined from P\ by identifying the tensor factor B s in the range of 
P 1 with B d ^' s C A d ^ x \ and 

(8) the total transition operator P : A c — > A c , defined by equation (2.3). 



To see that the operator in step (8) is well-defined, observe that the subcells 
are labelled by C x S, and that the sets R x = {(d(s; x) , s)\s G S} C C x S, which 
describe the range of P x are disjoint. Then, by Proposition IV.4.23(i) in [Tak], P 
is completely positive. Since PI = I, this implies that P is norm continuous, and 
consequently has a unique extension by continuity to the whole quasi-local algebra 
A c . 

Of course, when we think of a lattice, the injective maps d(s; •) will typically be 
lattice translations. If we choose all d(s; •) to be the identity, we obtain a system of 
non-interacting cells A x . Note that the subcell decomposition is only relevant in the 
range of P x , not in the domain. Thus in each step the subcell decomposition of the 
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previous step is obliterated. In particular, the second iterate P 2 of a QCA cannot 
be written in the same form: the algebras P 2 (A X ) do not commute with each other. 
This does not contradict the necessity of the commuting range condition explained 
above, because the product form of P 2 is also lost. Note that this is not an artefact of 
our quantum generalization: even in the classical case the second generation updates 
of a PCA are no longer independent, hence they no longer satisfy the definition of 
a PCA. 

The adjoint operator of P takes states into states, and we will usually denote 
its action by u> i— > uoP. It is easy to see (e.g., using the Markov-Kakutani Theorem, 
Theorem V.10.6 in [DS]) that any QCA has an invariant state, i.e., a state p such 
that p o P = p. A QCA is called ergodic, if there is only one invariant state for P. 
The main problem addressed in this paper is to find sufficient criteria for ergodicity 
in terms of the given local data P\ and d(s; •). We are also interested in stronger 
versions of ergodicity, e.g., the property that P n contracts in norm to the invariant 
state, i.e., 

lim \\P n (A) -p(A)l\\ = , 

for all A G A c . A further closely related problem is to estimate the decay of 
correlations in the invariant state. 

To see what is involved, it is good to look at the most trivial example: a non- 
interacting particle system. Then we have only one type of subcells, and d(l; •) = 
id£. P is simply the infinite tensor product of copies P x of a fixed operator Pi, acting 
on isomorphic finite dimensional algebras A x . It is obvious that the restriction of 
an invariant state for P to a single site is invariant for Pi, and, conversely, any 
product state formed out of invariant one-site states will be invariant for P (the 
latter construction need not be exhaustive). Hence P is ergodic if and only if Pi is 
ergodic. It is plausible that contractivity properties should also carry over. Assume, 
for example, that 

\\P?(A)-pi(A)l\\<e n \\A-pi(A)l\\ , 

for all A G A x , and p\ the unique invariant state of P\. Does this imply a similar 
bound for PI A direct estimate gives indeed a similar bound, but with e n multiplied 
by the number of sites (see Section 4.2). Hence this approach is not feasible on an 
infinite lattice. What one needs is a norm such that a contractivity estimate for a 
tensor product of completely positive operators P x is not worse than the maximum 
of the estimates for the factors. This is precisely the role of the oscillation norms 
used in the classical results of [LMS]. Their quantum analogue will be studied in 
the following two sections. We will then return to QCAs in Section 5. 
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3. Oscillation norms on C*-algebras 

3.1. Definition and basic properties 



In this section we want to generalize the notion of oscillation norm to the non- 
commutative setting. The basic idea remains the same: we consider some operations 
5 a , which annihilate constants, i.e., 5 a (T) = 0. In the classical case these operations 
measure the effect of the spin flips at different sites, i.e., 

(<W)(CT Q , Onota) = 2 (/( _(T a' °"nota) ~ /K^nota)) , (3.1) 

where a not a stands for all spin variables at sites other than a. Then we can say that 
an observable / is nearly constant, if the "oscillation norm" ^ a ||5 Q (/)|| is small for 
all a. 



For classical systems with more than two states per cell, as well as for any 
quantum system, there are many ways of "flipping" a single cell. There are several 
proposals in the literature how to take this into account. However, all proposals 
agree that the "total oscillation norm" of a lattice system should be the sum of the 
oscillations of each cell. Moreover, we will see below that the total oscillation norms 
in these different approaches are equivalent, whenever the local cells are described 
by finite dimensional algebras. Perhaps the simplest proposal [MZ] is to define the 
oscillation of an observable A localized in a single cell as 



\\\A\\\ = \\A- V (A)1\\ , (3.2) 

where rj is the normalized trace (or any other state) on the cell algebra A x . Other 
approaches use a family {5 a } of operators for each cell, and one can consider the 
"sup-oscillation norm" 

p| = sup 11^)11 • (3-3) 

a 

It is also suggestive to define 

Mh=sup mzmi , ,3.4) 

where the supremum is over all pairs of states, and d is some metric on the state 
space. In the classical case one could restrict the supremum to pure states, so that 
|||^4|||d is just the Lipshitz constant of A with respect to the metric d [Mae]. Finally, 
one may use for the single cell precisely the same form as for the total oscillation, 
namely a sum 

Pll|:=£ HM^)|| (3-5) 
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over "elementary" oscillations S a (A). This is the approach used by Matsui [Mai]. 
We will also adopt it, mainly because it agrees best with the subcell structure of 
cellular automata: the propagation maps will introduce a reshuffling of subcells 
between different main cells, and this process preserves oscillation norms only if the 
oscillation within each cell is defined by the same mechanism as the total oscillation. 
This approach also simplifies the presentation in the sense that we can use the same 
results about oscillation norms for the single cells as well as for the whole system. 



2 Definition. Let A be a C* -algebra, 1 G A. Let I be an index set and {5 a \a G /} 
a collection of bounded linear operators 8 a : A — > A, such that 

(1) 8 a (l) = 0, Vae/; 

(2) Af^ := {A G A\ Eae/ ||M^)|| < °°} ^ \\-\\-dense in A; 

(3) There exists a state 1] G A* (reference state), such that for A G Af^: 

\\A- V (A)1\\<J2\\UA)\\ ■ (3.6) 

Then, for A G A, 

P||| :=J2\\ 6 «( A )\\ 
is called the oscillation norm of A G A. 



Obviously, ||| • ||| is a seminorm on A. It satisfies 

inf||A-Al|| < \\A-r](A)l\\ < P||| . (3.7) 
Aec 

In particular, P||| = implies that A is a multiple of the identity. 

A simple argument shows that an infimum such as the one on the left hand side 
of (3.7) is attained at some A = X(A) in any normed space A with a fixed element 
1 6 A In a Hilbert space we can even assert that A ( A) is uniquely determined, 
and depends linearly on A. In a C*-algebra, however, the shape of the unit ball is 
different and although, for hermitian A, X(A) is uniquely determined, it is a non- 
linear functional. If the overall bound in (3.7) holds, however, we can easily make 
every state an admissible reference state, albeit for the oscillation norm P|||' = 
2p|||: if A is such that \\A- AI|| < P|||, then, for any state rj, \\A - rj(A)l\\ = 
|| (A - AI) + r]{A- AI)I|| < 2 || A - AI|| < 2p|||. This shows that the choice of the 
reference state is largely arbitrary. However, the existence of some such state is a 
non-trivial constraint on {5 a }, as the following example shows. 
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3 Example. Let A = ® fe z A 1 be the quasilocal algebra on the one-dimensional 
chain, with A 1 = A 1 for all i. Denote by r : A ^ A the translation automorphism, 
defined by 

T i® Ai ) =0 Al+1 ■ 

We define an operator 5 : A — > .4. by 

5(A):=r(A)-A . 

Then 5(A) = is equivalent to the translation invariance of A, which in the quasi- 
local algebra A is equivalent to A e CI. Thus the one-element collection {5} satisfies 
the first two conditions of Definition 2 (with Afm = A), and, moreover, 5(A) = 
implies A e CI. But condition 3 is violated in the strong form that there is no finite 
constant C such that (3.7) holds as \\A - AI|| < C|||A|||. 

To see this, pick some element Ax e A 1 \ CI in the one-site algebra, and a 
function / : Z -> H with /(i) > 0, and £»/(*) = L Then set 

i 

Let wi be a state on A 1 , and let w-j 2 = ® ie z u;i be the infinite product state on A. 
Then 

||A - AI|| > uif{A - AI) = wi(Ai) - A 

n^)ii = || X)(/co - /(< - < Piii £i/(o-/c-i)i • 

% i 

with a suitable choice of u\ and A we find \\A — AI|| > infy ||Ai — A'I||, and hence, 
with = \\8{A) || 



C" 1 < 



m < T\m-f(i-i) 



By choosing / to be slowly varying, e.g., f(i) = (1 — /z)/(l + /i) for |it — > 1, the 
right hand side can be made arbitrarily small. A 

Oscillation norms defined with a single 5 a are precisely those of the form (3.2), 
i.e., 5(A) = A — rj(A)l. The lower bound (3.6) is then equivalent to |||A.||| < |||A.|||. 
In the other direction we have the estimate (see Appendix IV of [MZ]) 

\\\A\\\ = J2\\^( A -V(A)1)\\< miUjPlllo • (3-8) 

of \ a / 

Hence, if there are finitely many oscillation operators 5 a , or the above sum is oth- 
erwise convergent, the norms ||| ■ ||| and ||| • |||o are equivalent. This is not a big 
surprise, since on a finite dimensional A, all semi-norms which vanish exactly on 
the constants are equivalent. 

In the next example we will consider a special case of Definition 2, which will 
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sometimes be especially convenient. The oscillation norms used in [Mal,Ma2, 
Ma3] and [MZ] are of this form. 



4 Example. We say that the operators 5 a satisfy Matsui's condition, when 

J2UA)=A-rj(A)l , (3.9) 

for some state r\. Then, by a simple application of the triangle inequality, condition 
3 of Definition 2 is satisfied. Another situation in which this condition comes up 
naturally is the following: Let A = Ai n be the algebra of n x n-matrices, and let G 
be a finite group. Consider an irreducible projective representation g i— > U g G M. n 
of G, and set 

6 g (A) = ±-(A-U g AU;) . (3.10) 

Then J2 g U g AU* commutes with all U g , and is hence of the form r}(A)l. 

Clearly, n is an invariant state with respect to all U g , and must therefore be the 
normalized trace of Ai n . Since for the identity element e G G we have 5 e = 0, 
it suffices to take the above 5 g for g G / = G \ {e}. With this choice, Matsui's 
condition (3.9) holds. Note that since ||A|| = the oscillation norm may be 

written in the suggestive form 

111^111 = iTJi Ell [^^ill • 

|G| 5 6G 

The simplest special case is to take the U g as the three Pauli matrices a a G M.2, 
and U e = 1: the product of any two of these operators is in the same set, up to 
a phase. The group G (\G\ = 4) consists of the rotations by n around the three 
Cartesian axes in 1R 3 . In this case one also finds easily that the reference state r\ is 
uniquely determined: let rj' be another reference state. Then 

1 3 1 
Iki - v'Ml\\ < III^IH = - J2 \\[vu<7i}\\ = 4(l|2^|| + ||2a 2 ||) = 1 

i=l 

implies that t]'(o~i) = 0. Repeating the same argument for the other components we 
find that n' = rj has to be the normalized trace. A 

The infimum on the left hand side of (3.7) is the standard quotient norm of 
A/&1. Hence ||| • ||| can be considered as a proper norm on the subspace ^Ifin/Cl of 
this quotient. We now show that this norm turns Afm/Gl into a Banach space. 

5 Lemma. /(CI is ||| • \\\-complete. 



Proof : Let {A n } C ^Ifin/Cl be a ||| • |||-Cauchy sequence. It follows that ||A n || < 
|||-An|||- Because A/ 'CI is complete for the quotient norm ||-||, there exists an A G 
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A/ CI such that \\A - A n \\ -> for n ->• oo. We have to show that A e An/CI. 
Let J' C J be a finite subset of the index set. Then, because each 5 a is bounded, 

V ||5 a (A)|| = lim V ||<5 a (A n 

ceEl' ael' 

< lim sup ||5 a (A T 

< lim |||A n |||. 

n — >oo 

This limit exists because A n is ||| • |||-Cauchy. Taking the supremum over all finite 
I' C /, we find that |||A||| < oo. Similarly, we find that 

V ||<y a (A-A n )|| < lim p m -A n ||| , 
aer 

and since this bound is independent of I' C /, we have that lim n \\\A — A n ||| = 0. 



Since we have not postulated any further properties of the operators 8 a , the 
space Afm does not come with a natural algebraic structure. However, in the special 
case, where S a (A) = i[A, D a ] are derivations, but also in the case (3.10), we get 
\\6 a (AB)\\ < \\5 a (A)\\ \\B\\ + \\A\\ \\5 a (B)\\, and hence 

\\\AB\\\ < \\\A\\\ \\B\\ + \\A\\ \\\B\\\ . (3.11) 

In particular, becomes a Banach algebra with the norm \\A\\ X = \\A\\ + A|||A|||, 
for any A > 0. 

The main reason for introducing oscillation norms is that in the case of large 
systems it is often easier to establish contractivity in this norm than in the C*-norm 
||-||. Nevertheless, as the following Proposition shows, this contractivity is sufficient 
to establish convergence of the iterates in the C*-norm, and hence ergodicity. Note 
that while the oscillation norms (3.2), (3.3), (3.4), and (3.5) are "equivalent" infinite 
dimensional situations, the best constant e in the assumption of the Proposition will 
depend on the choice of norm. 



6 Proposition. Let A be a C*-algebra with oscillation norm \\\ ■ |||, and consider a 
linear operator P : A — > A such that \\P(A)\\ < \\A\\, P(l) = J, and, for some fixed 
e < 1, and all A e Afm, 

\\\P(A)\\\ <e p||| . 
Then there exists a unique state p E A* , such that p o P = p. Moreover, 

(1) lim^oo \\P n (A) - p(A)l\\ = 0, for all A e A, and 

(2) \\P n (A) - p(A)l\\ < 2 e n P|||, for all A G Ag*- 
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Proof : Fix an element A e Af^ and consider the sequence A n = P n (A). Choose 
some numbers A n , for example A n = rj(A n ), such that 

\\A n ~ A n l|| < |||A n ||| 

Then {A n } is a Cauchy sequence: taking n > m without loss of generality, we get 

|A n — A m | < ||A n l — A n \\ + \\A n — A m l|| 

< |||A n ||| + \\P n ~ m || \\A m -\ m l\\ 

<(e n + e m )P||| . 

Let p(A) = lim 

m^oo A m . Inserting this definition into the previous inequality, we 

find 

\\ n -p{A)\<e n lA\\\ , 
and for A e A^ it follows that 

\\P n {A) - p{A)l\\ < \\P n {A) - A n I|| + |A n - p(A) 

< P n |||+e n |||A||| 

< 2e n p||| . 

p(A) is linear, because A i— > p(A)l is the limit of the linear operators P n . p(A) is 
also a ||-||-continuous functional on .4.^: 

\p(A)\<\\p(A)l-P n (A)\\ + \\P n (A)\\ 

<2e n \\\A\\\ + \\A\r^\\A\\. 
Let p denote the continuous extension of p on A. Then < ||A||, i.e., ||p|| = 1, 

and also = lim n P n (T) = I. Hence p is a state. 

For showing the convergence (1), let A e A, and 5 > 0. Then we may pick As G ^Ign 
with \\A - A s \\ < 5. It follows that 

\\P n (A) - p(A)l\\ < \\P n (A - A s )\\ + \\P n (A s ) - p(A s )l\\ + \p(A s - A)\ 
< ||^_^|| + 2e -p 5 ||| + ||^-^|| 

<25 + 2e n p 5 ||| . 
Hence, for sufficiently large n, the left hand side becomes arbitrarily small. 

Finally, suppose v G A* is another fixed point of the adjoint of P, i.e., v o P = v. 
Then 

v(A) = lim (z/ o P n ) (A) = u(p(A)l) = p{A) ■ u(l) . 
Hence, if v is also normalized, we have v = p. 



We want to use this criterion to establish ergodicity of cellular automata de- 
scribing spin systems. The operator P will then be built up from local operators P x . 
We thus have to construct an oscillation norm for composite systems, i.e., for tensor 
products of algebras with oscillation norm, and then have to apply Proposition 6 to 
the product system. 
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3.2. Tensor able oscillation norms 

We now take the algebra A = A c to be the quasi-local algebra over a lattice £, to 
each site of which is attached a unital C*-algebra A x ', i.e., 

A = A C = $$A X . 

xec 

Assume now that we are given an oscillation norm for the algebras A x at each site. 
We would like to assemble from this an oscillation norm for A c . 

The basic idea is very simple: as the collection of "flip" operators we simply 
take the union of the flip operators for each site, i.e., 

{J{5W\aeI*} , (3.12) 

where 6a = d x <8> id £ ^ x ^ is just the action of d x at site x of the quasilocal algebra. 
The oscillation norm of an element A e A is then 

iiw« = ££K x) wil • ( 3 - 13 ) 

This quantity is always defined, but possibly infinite. So we can define, as before, 

Aon = {A e A\\\\A\\\ < oo} . 

Clearly, the finite tensor products ® xe \A x , with A x e Af^, and their linear com- 
binations are in Af^. By Definition 2, Af^ C A x is norm dense, and by definition of 
the quasi-local algebra the finite tensor products span a dense subspace of A. Hence 
Afin fl A\oc, and a fortiori Abd., is dense in A. In order for ||| • ||| to become an oscil- 
lation norm, we need to establish the estimate (3.6) in the definition of oscillation 
norms, with the obvious candidate 

V = ®rj* (3.14) 

xec 

for a reference state. For this, we need the following "stabilized version" of (3.6). It 
is automatically satisfied in the commutative case (cf. Proposition 10 below). For 
the notion of complete boundedness, see the Appendix. When r\ is a state on the 
C*-algebra A, we denote by fj : A — > A the operator fj(A) = rj(A)l. 
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7 Definition. An oscillation norm defined by operators {5 Q |a G /} on a C*-algebra 
A, with reference state 77, is called tensorable, if each 5 a is completely bounded, 
and 

\\A-(id Mn ^)(A)\\<Y,\\(MM n ®Sa)(A)\\ , (3.15) 

ceEl 

for Ae M n (€) <g> A. 



At first sight it may seem rather special to allow tensoring only with M. n . 
However, as in the definition of complete positivity this is sufficient to give the 
analogous statement for all C*-algebras. 



8 Lemma. Let {5 Q |a G /} define a tensorable oscillation norm on a C*-algebra A, 
and let M. be any unital C*-algebra. Let M®Abe the minimal C*-tensor product, 
and A G M <g> A. Then 

||A-(id^®?))(A)||<^||(id M ®5 a )(A)|| . (3.15') 

ceEl 



Proof : We need the following basic observation: if p 1 is a net of Hilbert space 
operators with ||p 7 || < 1, converging strongly to the identity operator then, for any 
bounded operator A, 

limllpMpJI = \\A\\ . (3.16) 
7 " ' 11 

Indeed, the inequality limsup 7 ||p*Ap 7 || < limsup 7 ||p 7 || ||A|| < ||A|| is trivial. On 
the other hand, let 0, ip be unit vectors such that \(<p, Aip)\ > \\A\\ — e, and let 7 be 
sufficiently large such that \\p 7 4> — <f>\\ , Wp-yi^ — V'll ^ e - Then 

\\p;A Pl \\ > \(p^,Ap^)\ > \{<t>,A^)\ -2e\\A\\ > \\A\\ - e - 2e \\A\\ , 
and (3.16) follows by taking the inferior limit, and e — > 0. 

Without loss we may take M. and A to be faithfully represented on some Hilbert 
spaces. Then the minimal C*-tensor product M. <8> A is defined as the C*-algebra 
generated by operators of the form M <S> A, with M G M. and A G A, respectively. 
Let Pj denote a net of finite dimensional projections in the representation space of 
M, converging to the identity, and introduce the operators p 7 (M) = p 1 Mp 1 . Then, 
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for A e M <S> A, we have A = (p 7 <g> id^)(A) e -M n <g> A, where n is the dimension 
of p 1 . Hence, because (j? 7 ® id_4) and (id^ ®v) commute, 

\\{p^®id A ){A- (id M ®fi){A))\\ = \\A - (id>,„ ®rj){A)\\ 

<EIK id ^ ®^){A)\\ 

= Yl II fa ® id ^) ( id -^ 0( M 

<]T|K id ^®^K A )ll • 

Hence the result follows, by applying (3.16) to the left hand side of this inequality, 
and the family of projections p 1 <S> I. 



The basic result concerning the tensor product of algebras with oscillation norm 
is the following. The second part was proven in a special case by Matsui [Ma4]. 



9 Proposition. Let {A x \x E £} be algebras with tensorable oscillation norms, 
defined by {S x \a e an d with reference states rf e (A x )* . Then 

(1) A = <S) xe cA x JS an a lg e bra with tensorable oscillation norm given by (3.12), 
with reference state (3.14). 

(2) Aftn fl A\ oc is dense in Af^ with respect to ||| • |||. 



Proof : We show (2) first. For any finite subset Ac£, define 

x<£A 

where f} {x \A) = [f) x ® id A{x} ) (A) . Therefore the range of fj A is contained in the 
local algebra <S> xe \A x . For A e A\ OC: we trivially have 

\\fi AC (A) - A\\ A -^0 , (*) 

because A absorbs the localization region of A. Consequently, because A\ oc C A is 
|| -||-dense, (*) is valid for all A e A. 

We show next that 

\\\^(a)-a\\\ a ^o , 
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for A G Asn. One has 

^ c (a)-a\w=y: j:\k x) (v ac -^)(a)\\ 

xEA ael x 



E E IIW-a £ )wil + E E llW-u^wil 



x£A a6i" a 



A c 



= E E iitf 

because for x G A the operators r) A 
on the factor l (x) , so 5i x) o r? AC = 0. 



x^A ceEl x 

id £ )^(A)|| + ^^||^)(A)|| 

x^A ael x 

and 5a commute, and for x A, 5a" y acts 
In the limit A — > £ the second term vanishes 



because the sum over all x G £ converges for A G ^4fin- The first sum is termwise 
dominated by the sum defining 2|||A|||. However, because 6a (A) G A, each term in 



this sum goes to zero by (*). Hence the sum goes to zero by dominated convergence. 



Hence \\\rj A (A) - A\\\ -> 0, and Afm H A\ c 



is 



-dense in Aon- 



To prove (1), we have to verify Definition 2. The boundedness of the 8 a x ^ follows 
from the complete boundedness of the operators S x on A x . The property 5^(1) = 
is evident, and we argued for the ||-||-density of Af^ above, before (3.14). Hence only 
the estimate (3.6) remains to be seen. We begin by showing it for A G A^ H Ai oc , 
say A localized in a finite region Ac£. We conveniently label the sites in A as 
1, . . . , |A|. Then with a telescoping sum we find 



V X (A)\\< 



xeA 



|A| 

E 

k=l 

|A| 



fj x <g> (id fc -fj k ) 



x<k 



<g) id y (A) 

y>k 



< 



E 

k=i 

EE 



x<k 



\A-f)W(A)\\ 



keA a ei k 

where i)^ : A ^ A is the operator ff ® id £ ^' . At the last estimate we used that 
the product of the rf x > is a completely positive unit preserving operator, which has 
hence norm 1, and, of course, the tensorability of ||| • |||. Hence the required estimate 
holds for A G Asn fl A\ oc - In particular, it holds for the approximants r) AC (A) of a 
general A G Afm- Since by part (2) of the Proposition these elements approximate 
A in both norms, both sides of the estimate converge as A — > C This completes the 
proof that ||| • ||| is an oscillation norm. It is tensorable, because we can include an 
additional tensor factor M. n in the product defining A, and use the same estimates 
as above to establish (3.15). 



With this Proposition the construction of oscillation norms for infinite systems 
is reduced to the construction of tensorable oscillation norms for the one-site alge- 
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bras. Recall from (3.8) that all oscillation norms for which J2 a ei W^aW < 00 are 
equivalent. By summing this estimate over all cells we conclude similarly that all 
oscillation norms on a composite system, for which 

E ra cb <c<oo (3.17) 
aei x 

with a constant independent of x, are also equivalent. In particular, this holds 
for the norm in [MZ]. Equivalence can also be shown for global oscillation norms, 
which are defined as the sum of local oscillations, which are in turn given by some 
supremum (see (3.3) or (3.4)) [Mae]. 

Showing tensorability is especially easy in the classical case: 



10 Proposition. On an abelian C*-algebra A, every oscillation norm is tensorable. 



Proof : We can set A = C(X), for some compact space X. Let M be any C*- 
algebra. Then M. ® A = C{X,M), the algebra of .M-valued continuous functions 
on X , with the norm 

H|:=supP(x)|| = sup|$(A(x))| , (3.18) 

X x,$ 

where the supremum over $ is with respect to all linear functionals on M. with 
|| $|| < 1. Then ((id^ <g)5 a )^4)(:r) = f S a (x,dy) A(y), where we have used 5 a to 
denote both the operator and its integral kernel. Applying the estimate (3.6) to the 
continuous function x 1— > $(A(x)), we obtain 

\\A-(id M ®fj)(A)\\ = sup 

x,<& 



< 



supE 



*(A(x)) - J V (dy)*(A(y)) 
8 a (x, dy) *(A(y)) 



sup 



ceEl 



< E SU P SU P / 5 a (x,dy) $(A(y)) 
a ei * x J 

= Y,\\QdM®8a)(A)\\ . 



In the non-commutative case, tensorability is a non-trivial constraint on oscil- 
lation norms. The following is a handy criterion. 
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11 Lemma. Suppose that the operators {5 a \a £ 7} defining an oscillation norm 
HI • I on a C*-algebra A satisfy Matsui's condition (see Example 4 above). Then it 
is tensorable. 



Proof : Let A e M n (<D) <g> A. Then, by (3.9), 

^(idx„ ®6 a )(A) = (id Mn ®J2^)(A) = (id Mn ®(id A -fi)yA) 

= A - (id Mn ®f^(A) . 
By taking norms on both sides we find (3.15). 



Finally, we record for later use that for tensor product operators a version of 
(3.11) holds without further assumptions: using that in a minimal C*-algebra tensor 
product ||A(8)5|| = ||A|| ||B||, we find 

\\\A® B\\\ = \\\A\\\ ||S|| + \\A\\ \\\B\\\ . (3.19) 
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4. Contractivity of tensor product operators 

4.1. Oscillation norm estimate 

In the simplest, non-interacting case of a cellular automaton the total transition 
operator P is the infinite tensor product of the one-site transition operators P x . If 
we know that each P x contracts exponentially with rate e < 1 to a multiple of the 
identity, can we also assert this about P? This turns out to be the crucial question 
for developing ergodicity estimates for quantum cellular automata. We will show in 
this subsection that, provided we use oscillation norms to express contractivity, the 
product P indeed contracts with the same rate. In fact the validity of this bound 
is the main reason for considering oscillation norms. In order to make this point 
more precise we show in the next subsection how other norms fail to give the desired 
estimate. 

The most straightforward definition of contractivity of a transition operator P 
is given by the estimate |||P(A)||| < e |||A||| as in Proposition 6. The best constant e in 
this estimate is conveniently denoted by |||P|||; it is the norm of P as an operator on 
Aon. However, this quantity is not appropriate for the study of composite systems, 
since even in the classical case we may have ||| \&m ®P\\\ > \\\P\\\ (cf. the example at 
the end of this subsection). Therefore, we use a version of |||P|||, which is "stabilized" 
with respect to coupling the system to an outside world. 



12 Definition. Let A and B be algebras with tensorable oscillation norms gene- 
rated by {S a \a E 7} and {5p\/3 £ J}, respectively. Let P : A — > B be linear and 
completely bounded. Then the completely bounded oscillation norm of P, 

denoted by |||P|| c b> JS defined as the smallest constant e for which the inequality 

I OdM„ ®&P))(A)|| < 6 II (^M n ®S a )(A)\\ (4.1) 
/3eJ aei 
holds, for all n e IN, and A e M n ® A. 



In particular, for n = 1, we get |||PA||| < |||P||| c b |||A|||, i.e., |||P||| < |||P||| c b- 
Precisely as in Lemma 8 one sees that if the bound of the form (4.1) holds for all 
algebras M. n ®A, it also holds for all minimal C*-tensor products M. ®A with other 
C*-algebras. The crucial property of this norm is given in the following Theorem. 
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13 Theorem. Let P x : A x — > B x , x G C be a family of completely positive unit 
preserving operators between algebras with tensorable oscillation norms. Let A = 
<S) xeC A x and B = <S> xeC B x be equipped with the oscillation norm described by 
Proposition 9, and let P : A — > B be defined by P = (g) xeC P x . Then 

\\\P\\\cb = sup |||P x |||cb • (4-2) 



Proof : The inequality |||P|| c b > su Px |||^a;|||cb is trivial, because the estimate (4.1), 
written for P, and an observable A G 1m ® A x C 1m ® A reduces to the corre- 
sponding estimate for P x . 

For the opposite inequality we have to show that, provided each P x satisfies the 
estimate (4.1) with the same constant e, i.e., e > |||P x |||cb for all x G £, then so does 
P. Consider A G M <E> A, and one term in the sum on the left hand side of (4.1), 
say for x G £, and f3 G J x . We have to estimate 

(icU, ®{^ ] P)) (A) = (id M ® (g) P y <8> id^) (id M <g> (g) id^ H ®(5%P x j) (A) . 

y=£x y^x 

The first parenthesis is a contraction because each P y is completely positive and 
unital. Hence taking norms, summing over (3 G J x , and using Definition 12, we find 

£ lOdx^^p))^)) 



|PJ £ (id M ®(g)id^^)(A) 



<e£ II (id^®^))^) 



The sum of these inequalities over all lattice points x G C is the desired estimate, 
showing |||P|||cb < €. 



Hence, for a non-interacting system, it suffices to show |||P x |||cb < e < 1 to 
conclude ergodicity from Proposition 6. However, the explicit estimate of |||P x |||cb 
may still be a difficult problem. One way to handle it is the following decomposition 
property. 
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14 Lemma. In the setting of Definition 12, suppose that for each (3 G J, we have 
a decomposition 

with G(3 a : A — > B linear and completely bounded, and the sum strongly convergent 
on A. Then 

||| j P|||cb<sup^||G' /3Q || cb , (4.3) 
where ||-|| cb is the completely bounded norm. 



Proof : The proof is obvious by inserting the decompositions of SpP, and using the 
triangle inequality. 



If A and B are finite dimensional, such decomposing maps Gp a always exist. 
In fact, the necessary and sufficient condition for an operator (here 5pP) to allow 
a decomposition with given 5 a is that ker^P D f] a ker5 a . However, the right 
hand side of this inclusion is equal to CI by part 1 of Definition 2, and, clearly, 
CI C kerdpP. We conclude that, on finite dimensional algebras with oscillation 
norms defined by finitely many S a , we have |||-P||| c b < oo for all transition operators. 

If |||P|||cb = 0, we must have 6p(P(A)) = for all 0, so that P(A) = P U (A) = 
u)(A)l for some state uj on A. Near maps of this form we find transition operators 
with small oscillation norms: these could be called transition operators with large 
disorder, since in one step they wipe out nearly all memory of previous states. It 
is straightforward to see from the definition of the completely bounded oscillation 
norm that, for < A < 1, the equation 

|||(1 - A)P W + AP|||cb = A|||P||| cb (4.4) 

holds. (In fact, the inequality "<" already follows from the convexity of ||| ■ ||| c b)- 
Hence, for sufficiently small A, the ergodicity criterion Proposition 6 applies, and 
(1 — A)P^ + A.P has a unique invariant state, which will be close to, but not equal 
to UJ. 

In some cases, one can use the freedom of adapting the operators 5 a to the 
problem at hand to give a simple estimate of |||P||| c b- An example is the following: 
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15 Lemma. Let P : A — > A be a transition operator on a finite dimensional C*- 
algebra, and suppose that P is diagonalizable, i.e., it has a representation in the 
form 



N 



P = D + J2 X « D * 



a=l 



with D a Dp = SapDu, for < a, (3 < N , and Dq(A) = u(A)l for some state u> on 
A. Define an oscillation norm by setting 5 a = D a for a = 1, . . . , N. Then 

\\\P\\\cb = max ||A Q | a = 1, . . . , ivj . 



Proof : Note that the oscillation norm so defined satisfies Matsui's condition (3.9), 
because J2^ =Q P J a = id, and is hence tensorable. The lower bound |||-P||| c b > 
max|Ao,| follows by inserting the eigenvectors into the estimate defining ||| • ||| c b, 
and the upper bound follows from Lemma 14 with G^ Q = A Q 5 a/ 3id. 



In the classical case we may simplify the definition of ||| ■ ||| c b, by considering 
couplings to classical systems only. For Ising systems, i.e., the case considered in 
[LMS] and other principal papers on the subject, the stabilization can even be 
omitted entirely. 



16 Lemma. Let A, B be abelian algebras with oscillation norm, and P : A — > B a 
completely bounded linear map. Then 

(1) |||-P|||cb is the best constant e such that the estimate 

\\(id M ®(hP))(A)\\ < e W( [dM 

/3eJ a€l 

holds for all abelian C*-algebras M. 

(2) if B is two-dimensional, or the oscillation norm on B is defined by a single 5, i.e., 
J has only one element, then |||P||| c b = |||-P|||> J '- e -j the best constant is already 
achieved by taking M. one-dimensional. 



Proof : As in the proof of Proposition 10 the quantum observable algebra is reduced 
to a classical one by evaluating in appropriate states. Let A = C(X), B = C(Y), 
and let M. denote the C*-algebra of bounded functions on the index set J. Suppose 
the bound in (1) holds with this abelian algebra M, and let A e M. n <g> A. Thus, 
for each (3 e J, (id.M„ ®8pP){A) is a continuous .M n -valued function on Y. By the 
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definition (3.18) of the norm in M n <S> B there is a linear functional on Ai n , of 
norm < 1, such that 

(id Mn ®6pP)(A) = (^ (3 ®~S (3 P)(A) = 5 p P($f3®id A )(A) 

M n ®B B B 

where we used that $^ : M. n — > C, and C ® B = B. We now introduce the function 

A e M ®A, defined as A(J3,x) = §p(A(x)), or A(J3, •) = ($3 <g> id A )(A). Then the 
norm on the right hand side in the above equation is smaller than 

6 /3 P($ l3 >®id A )(A) 



sup 



B 



(id M <S>S^P)(A) 



< 



M n ®B 



J2 \\( [d M ®5 a )(A) 



Summing over /3, and applying the given inequality for Ai, we find 
J2\\(id Mn ®SpP)(A) 

and the result (1) follows, because 
(id M ®5 a ){A) 



sup II ($/? 



< sup 



*8 a ){A)\\ A 
^)(A)\\ A = 



id 



Mr 



>8a){A)l 



To prove (2), note that, on a two-dimensional algebra £>, all operators with 
<5/3 (1) = are proportional, so we may replace the definition of the oscillation norm 
by an equivalent one with \J\ = 1. Hence the claim follows from the observation 
that we used (1) only with the now one-dimensional algebra M. = C(J). 



However, in classical systems with more than two spin values per site we may 
have strict inequality |||P||| c b > I II P I II- I 11 the following Example 17 we even have 
|||P|||cb > 1 > |||-P|||- Hence P contracts exponentially to its fixed point, and the same 
is true for a non-interacting QCA with this one-site transition operator. However, 
for a system with non-trivial propagation maps this information is not sufficient, 
and only |||P||| c b gives a bound which is independent of the propagation. 

For computing |||P||| c b Lemma 14 may be helpful even in the classical case, since 
then the norms ||G/3 a || cb can be replaced by ordinary norms. 

We remark that with some of the modified definition of the oscillation norms 
on the single site observable algebras described at the beginning of Section 3.1, the 
stabilization can be avoided altogether in the classical case [Mae]. However, as 
already remarked in that context, this would be in conflict with our technique for 
proving the ergodicity criterion in the general quantum case. 
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17 Example. We take A as a system of two Ising spins, with its standard oscillation 
norm (3.1). This norm can be written as 

1 max |/(-0i, <j 2 ) - /(<Ti, cr 2 )| + J max |/((7i, -cr 2 ) - /(<7i, 2 )| 



0"1,0"2 



2 0"1,0"2 



fl,CT2 



X 0-10- 2 /(<Ti, (T 2 ) 



0"1,0"2 



where the variables 0^ are ±1. From the second form it is obvious that the unit ball 
of I ■ I is the cartesian product of 1R (corresponding to multiples of the identity), 
and an octahedron in 1R 3 . We now consider a transition operator P : A — > A, given 
by the matrix 

/0 1 5 6 ^ 
2 10 
7 5 
\8 4 / 
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(4.5) 



in a basis in which the components of / are (/(++), /(H — ), /( — h), /( ))• The 

oscillation norm |||P||| is readily computed, by applying P to the 6 extreme points 
of the octahedral unit sphere of ||| • |||, and computing the oscillation norms of the 
images. The result is 

iii^iii = § • (4-6) 

We now couple the system to an additional Ising spin, denoted by 00, and described 
in the algebra M. = C . What we have to estimate is the operator norm of 
(idx ®P) : M ® A — > M ® A with respect to the norm 

\\\f\\\M®A-=\ max 1/(00, -01,02) - /(0O, 01, 02)| 

Z O"0, cr lj°"2 

+ \ max 1/(00,01,-02) - /(0o, 01, 02)| • 

This norm is not characterized as easily as before. The unit ball has two unbounded 
directions, and the compact convex set in the remaining 6 dimensions is bounded 
by 48 hyperplanes. We did not succeed in computing all the extreme points of 
this polytope, so we have no explicit expression for ||| idx <S>-P|||. However, any 
expression \\\{id M <3P)f\\\ M ® A with |||/|||m®^ = 1, is a lower bound on |||P||| c b- 
Taking for / one of the 28 extreme points of the unit ball known to us, namely 

/ = (/(+ + +),-..,/( )) = (2, 1, 1, 0, 0, -1, -1, 0), we find 

13 

|||P|||cb> III id.M <g>P||| > — • (4.7) 

A 
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4.2. Estimates in other norms 



In this section, which is not needed later in this paper, we show how an estimate of 
the form (4.2) fails, if we use some criteria different from oscillation norms to define 
contractivity of the factors. To us this is the main motivation for using oscillation 
norms in the first place. 

Since P x l = I, "contractivity" has to be defined in a way ignoring this known 
fixed point. A natural approach is to consider contractivity in the quotient norm 
of A/Gl, i.e. 

\\A\\':= mf |L4-AI|| . (4.8) 
Similarly, for P : A — > A we define 

||P||' := sup{||PA||' ||A||'<l} , (4.9) 

and call P "contractive", if \\P\\' < 1. The following example shows what kind of 
estimate we can expect for this norm of a tensor product of transition operators. 



18 Example. We consider finite classical systems (A = C(O), O a finite set), for 
which any transition operator is of the form (2.2) 

P/(u/) = X>(u/,»7)/fa) • 
v 

One easily checks that 

||P||' = imaxV|p(a;,77)-p(u/,77)| . (4.10) 

Z u),oj' — 
V 

From this formula, and the positivity and normalization conditions for the tensor 
factors, one easily finds an estimate for tensor products, namely 

<g)p x \\ <J2w p *W ■ ( 4 - n ) 

xeA xeA 

Note that this estimate grows with the size of the region A, and becomes completely 
useless for infinite regions. Hence the question is whether this trivial estimate can 
be improved upon. 

As a simple counterexample, consider an Ising spin system (i.e., Q x = {+, — } 
at each site x), and all factors P x = Pi equal. Let x+ an d X- denote the functions 
which are 1 on the points "+" and "— ", respectively, and zero otherwise. Pi is 
characterized by the two probabilities 

p ± =p(±,+)=(P LX+ )(±) ■ 
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From (4.10), \\Pi\\' = \p+ — p-\. For the products x+ = (SWa X+ and P A = 
<S> xeA P x over a finite set A of N sites we get ||x+||' = 1/2, and 

\\p A xi\\'>l\(p A xi)(+---+)-(p A xi)( ) 

= ^\p+ -p"\ > ^\p+ -p- \ iVmin^.,^}^" 1 

= A^min^-,^}^" 1 ||Pi||' \\x+\\' ■ 
Picking both p + and p- close to 1, we see that ||-P A || / niay come arbitrarily close 
toAr||Pi||'. A 

Of course, the bound (4.11) cannot be improved upon in quantum systems 
either, and the norm (4.10) remains useless in infinite quantum systems, as well. 

Another alternative to oscillation norms, which seems plausible at first sight, 
is to use the observation that for norms of Hilbert space contractions with a known 
fixed point the right hand side of the analogue of (4.11) can be improved to a 
supremum. This suggests the use of the Hilbert space norms 

\\A\l = {u{A*A)-\u{A)\ 2 ) 1/2 (4.12) 

on A, and its associated operator norms for some state to. It follows from the 
complete positivity of transition operators that P is a contraction with respect to 
this norm, provided that uj is invariant under P (i.e., uj o P = uj). This may not 
seem like a severe restriction, since we know that any transition operator admits 
an invariant state. We then define ||-P||^ as the best constant in the inequality 
\\P(A)\\l <e\\A\\' u . The inequality 

holds, and it seems that we achieved our goal of finding a quantity that behaves well 
under composition. However, there are several drawbacks. First of all, the invariant 
states uii have to be explicitly known in order to compute any norm. Secondly, the 
ergodicity statement one gets from the inequality ||P|| W < 1 is rather weak: it allows 
no conclusion about states of the infinite system which are singular with respect to 
uj. Perhaps the most severe restriction, however, is that the propagation operators 
introducing interaction into the QCA setting by mixing different cells also fail to be 
contractions with respect to this norm, so the approach based on (4.12) seems to be 
limited to the trivial, non-interacting case. 

A similar criticism applies to the idea to use the spectral radius of P as 
an operator on A/ CI, denoted by cr'(P). Again, we have equality a' (Pi <g> P2) = 
maxj a' (Pi), but we have no control over this quantity for a product of two transition 
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operators, which we need to introduce interaction (see the proof of Theorem 19). 
As an elementary example consider, as in Example 17, a classical system with two 
subcell types of one Ising spin each, with the transition operator 



Pi = 





i 





° \ 





i 











i 








Vo 





1/2 


1/2/ 



In a system consisting of two cells we choose the propagation map d(l; •) for the 
first subcell type to be the identity, and the map d(2; •) for the second subcell type 
to be the flip. Then P\ contracts exponentially with rate 1/2 to its invariant state, 
which is the pure state on the configuration (H — ). However, the total transition 
operator P has three invariant states. 



Another conceivable alternative are the "sup-oscillation norms" introduced 
in (3.3), taken now not only as a way to define the oscillation norm in the subcells, 
but as a principle to construct the total oscillation norm. Defining the operator norm 
|||-||| cb in analogy to Definition 12, it is easy to show the analogue of Theorem 13. 
However, this does not suffice to give an_ergodicity criterion, since the estimate (3.6) 
fails for such norms, and consequently |||-|||-Cauchy sequences need not converge in 
A. A simple example demonstrating these claims is the sequence of averages of 
Ising spins over an increasing sequence of regions. The sup-oscillation norm of such 
averages goes to zero like the inverse number of sites in the average, but, of course, 
the sequence of averages is not convergent in the quasi-local algebra. 



5. Applications to Cellular Automata 

5.1. Ergodicity 

In order to apply the results of the previous section to interacting QCAs, we need 
tensorable oscillation norms on the algebras B s belonging to each subcell type using, 
say operators 5^ : B s — > £> s , a G I s . Then by Proposition 9 we have tensorable 
oscillation norms on each A x , and, consequently, on the algebra A c of the whole 
system. By Proposition 6 ergodicity follows from the estimate |||P||| < 1 for the total 
transition operator P. Thus we arrive at the following criterion. 

19 Theorem. Let a quantum cellular automaton be given according to Definition 1, 
and suppose that each B s is equipped with a tensorable oscillation norm. Then 

lllcb = |||-Pl|||cb • 

Consequently, if |||Pi||| c b < 1> the QCA is ergodic, i.e., there is a unique P-invariant 
state p on A c . 
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Proof : Let P = <S> xeC Pi be the infinite tensor product of the operators Pi acting 
in each A x separately. This is also the total transition operator of the QCA with 
the same Pi, but each d(s; •) equal to the identity. Consider also the automorphism 
D : A c — > A c which takes £ (x ' s) into B d(s ' x) ' s . Then 

P = D P . 

According to Theorem 13, |||P||| c b = |||Pi|||cb- Moreover, since the oscillation norm 
on B d ( s '> x )> s i s defined by the same operators 5^ as in B x,s , D is a ||| • |||-isometry, and 

|||D|||cb = !• Hence |||P||| c b = |||-D|||cb III^IIU = |||-Pi|||cb- 



Note that similarly to the non-interacting case the criterion |||Pi||| c b < 1 is a 
condition of "large disorder", which is satisfied as soon as Pi is sufficiently close 
to a map of the form Pi (A) = u>(A)l (compare (4.4)). A remarkable feature of 
this criterion is that it does not depend on the propagation maps d(s;-), which 
distinguish an interacting QCA from a non-interacting one. 

One might expect from Lemma 15 that, with a suitable choice of oscillation 
norms, |||Pi||| c b can be made equal to the largest modulus of eigenvalues of Pi apart 
from 1. However, this is not the case, since it is crucial for the proof of Theorem 19 
that each oscillation operator 5^ acts in only one subcell B s , and so the propagation 
automorphism D becomes a ||| • |||-isometry. Clearly, this property cannot be expected 
of the eigenprojections of Pi. Still, the second largest modulus of eigenvalues of Pi 
is always a lower bound to |||Pi||| c b- 
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5.2. The classical case 

Since in the usual definition of PCAs no subcell decomposition is used, it is not 
obvious how the classical results [LMS,MS1] can be subsumed under Theorem 19. 
In this section we show how this can be done, pointing at the same time to a possible 
generalization of Definition 1. 

The basic operator Pi defining the one-site transition probability of a PCA (see 
(2.2)) maps the one-site algebra A into the tensor product A = A® n , where n is 
the number of cells influencing the state in a single cell of the second generation, n 
is often finite, but we do not need this fact. We will construct a QCA, whose one- 
site algebra is A, with all subcell algebras B s isomorphic to A. Note that all these 
algebras are now commutative, so the "Q" in QCA only refers to the fulfillment of 
Definition 1. 

The total PCA transition operator P (see (2.3)) can be decomposed into three 
factors: the first is simply the infinite tensor product Pf of the operator Pi, mapping 
A c to A c . The information about the cell y to which a subcell (x, s) in the latter 
algebra belongs in the next time step is encoded in propagation maps d(s; •):£—>£ 
as before, i.e., y = d(s;x). This defines an automorphism D of A c as in the proof 
of Theorem 19. The final step is the sitewise application of the multiplication map 
M : A -> A, defined by 

n 

M(/i <g> • • • <g> / n ) = JI £ 

i=i 

Hence we get the factorization 

P = M £ D Pf , 

where M r denotes the tensor product of the copies of M acting at each site. 

It is precisely the use of the multiplication map M, that is the specifically 
classical element in this construction. M is also called the n th order diagonal of the 
algebra A because, writing A = C(O), and A® n = C(O n ), as we may for an abelian 
algebra, we have 

Mf(x) = f(x,x,...,x) . (5.3) 

Clearly, M is a *-homomorphism, and M(I<E>---<E>/<8>---<8>!) = /, for all positions 
of the factor / in the tensor product. The existence of the diagonal characterizes 
abelian algebras: if a homomorphism M of this description exists in a C*-algebra A, 
the commutativity of the tensor multiplication implies the commutativity of multi- 
plication. The adjoint of the diagonal map is a "state duplication map" , producing, 
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(5.1) 



(5.2) 



from a state on A, n copies of A in the same state. Its non-existence in the quan- 
tum case is the basis of "quantum cryptography" (see [EHPP] and references cited 
there). Of course, we can formally define M by (5.1), even in the non-commutative 
case. However, if A = Aid, n = 2, and $ is the unitary permutation operator 
exchanging the two factors, ||M(<E>)|| = n, i.e., ||M|| > n. Hence, on an infinite 
dimensional algebra, M is typically unbounded. But even in the finite dimensional 
case, ||M|| > 1 makes the definition of the infinite tensor product in (5.2) 
impossible. 

It is now easy to modify (5.2) so that we get a QCA in the sense of Definition 1. 
Its transition operator is 

P := D Pf M c . (5.4) 
One easily verifies that this is the QCA with one-site transition operator 

Pi : A -> A 
Pi = PiM . 

With J : A -> A, defined as Jf = f <g> l®^" 1 ), we have MJ = id A , and hence 

P N = M c P N J c , (5.5) 

for every power > 0, i.e., the PC A can be recovered completely from the QCA 
picture. 

For estimating the contraction rates of these operators we need the following 
Lemma. 



20 Lemma. Let A = C(O) be a finite dimensional abelian C*-algebra with an 
oscillation norm defined by operators 5 a of the form 

8 a f(<r) = c a (f(a a (a)) - f(a)^j , 

where c a G M, and a Q : O — > O. Then, for every n, the n th order diagonal M 
satisfies the estimate |||M||| c b < 1. 



Proof : Aid ® A can be identified with the algebra of .Md-valued functions on 
n, equipped with the norm \\g\\ = sup^ \\g(a)\\. Then, for / G Aid <8> A® n and 
a' = a a (a), the expression 

Pfc = f(<r', °, -/(o- 7 , ...,<?', a,..., a ) 



k times n — k times k — 1 times n—k+1 times 
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is bounded in the norm of Ma by jc,^ 1 1 1| (id>t d ®§a )/||- Hence 

|| ((id Md ®6 a M)f) (a) || = \c a \ \\f(<j', . . . , a') - /(a, . . . , a) \\ 

n-l 

= |c Q | I 

fc=l 

n-l 

<E||( id ^^ fe) )/| • 
fc=i 

Hence, after taking the supremum over a, summing over a, and using the definition 
(3.13) of the canonical oscillation norm on A® n , we get (4.1) with e = 1. The 
constant cannot be better than 1, because MJ = id^, and, obviously ||| J/||| = |||/||| 
for all /. 



Hence, from the factorization (5.5) we get |||-P||| c b < |||M ||| c b|||-P|||cb|||<^ III cb < 
|||P||| c b. From (5.4) and because D is an oscillation norm isometry, |||P||| c b ^ lll-f] C |||cb- 
Because P-f = D~ 1 PJ C , the last inequality is actually an equality. By Theorem 13, 
III-Pl |||cb = Ill-Pilllcb- Summing up these estimates, we have 

|||P|||cb < Wlcb = |||Pl|||cb • 

This is exactly the bound given in [LMS]. When comparing these results, however, 
note that Proposition 6 gives [^^(A) — p(A)l|| < 2e n |||A|||, without the superfluous 
factor (1 — e) _1 , which is present in [LMS]. 

Finally, we wish to point out that (5.4) also points to a possible generalization 
of the notion of QCA, which does not use subcell decompositions: we only have 
to replace M by some completely positive unital operator from A® n — > A. For 
such systems our method for obtaining oscillation norm estimates would apply un- 
changed, but they would no longer satisfy the condition of commuting ranges (see 
Section 2). In this sense the cells would no longer be "independently updated". 
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5.3. Decay of correlations in the invariant state 



We assume now that the ergodicity criterion |||Pi||| c b < 1 holds. What can be said 
about the unique invariant state p to which P contracts? It is clear that in the non- 
interacting case (i.e., d(s;-) = id), but also if Pi (A) — u>(A)l (i.e., |||Pi||| c b = 0), p 
will be a product state. Therefore, it is reasonable to expect that if |||Pi||| c b is small, 
we should obtain a state with good clustering properties. Moreover, in contrast to 
the ergodicity criterion, the propagation maps d(s; •) should enter the estimate for 
the correlation functions. 

We will first describe the relevant geometric properties of the d(s; •). For A C £, 
we will set 

ds(A-) = {d(s;x)\s G S, x G A} C C . 

This is the set of cells to which the interaction can spread from some site in A in 
one step. Similarly, we define dg(A) as the n th iterate of ds- For Ai, A 2 C C we 
define the correlation distance as 

c(Ai, A 2 ) = max {n G N|rfg(Ai) n d s (A 2 ) = 0} , 

i.e., as the last time step under which the two regions remain independent. When 
the maps d(s; •) are translations, it is clear that for large separation parameters r, 
the correlation distance c(Ai, A 2 + re) will asymptotically be proportional to r, but 
with a constant depending on the direction e. In this sense the following Proposition 
gives exponential clustering with a rate depending both on the direction, and on 

|||-Pl|||cb- 



21 Proposition. Let a quantum cellular automaton be given according to Defini- 
tion 1, and suppose that, with respect to some choice of tensor able oscillation norms 
on each B s , the ergodicity criterion |||Pi||| c b < 1 is satisfied. Let p denote the unique 
state such that p o P = p. Let 

A 1 g A Al = (g) A\ A 2 g A A2 = (g) A x , 

where Ai fl A 2 = are disjoint finite subsets of C. Then 

(1) P k (A 1 ) C ^4(Ai) f or a u k G N, 

(2) P k {A l ®A 2 ) = P k {A 1 )®P k {A 2 ), for all k < c(A l5 A 2 ). 

(3) |p(A 1 ®A 2 )-p(A 1 )p(A 2 )|<2(|||P 1 ||| cb ) c(Al ' A2) (pi||| P 2 || + Pi|| p 2 |||). 
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Proof : The first two statements are obvious from Definition 1 for k = 1, and follow 
for other k by induction. Then using Proposition 6 and Theorem 19, we get, for all 
k < c(Ai, A2), the estimate 
\p(A 1 ® A 2 ) - p(A l )p(A 2 )\ = \\p(A 1 ® A 2 )l - p(A 1 )p(A 2 )l\\ 

< \\ P {A 1 <g> A 2 )l - P k (A 1 ® A 2 )|| + ||P fc (^i ® A 2 ) - P fc (A0 ® P k (A 2 ) \\ 
+ \\P k (A 1 )-p(A 1 )l\\ 11^^(^2)11 + ^(^)1 \\P k (A 2 )-p(A 2 )l\\ 

< (|||Pi|||cb) fc Pi®^2||| + (|||Pi|||cb) fc P2|| lll^lll + (|||P 1 ||| cb ) fc ||A 1 || p 2 |||, 
where at the last step we used that ||P fc (A 2 )|| < ||A 2 ||. The result then follows from 
equation (3.19). 



Appendix: Complete boundedness 

To motivate the necessity of considering complete positivity and complete bounded- 
ness of operators on non-commutative C*-algebras we consider a standard example: 
the operator P : M n — > M n of transposition on the algebra of n x n- matrices. This 
preserves positivity and the identity element, and therefore seems to be a candi- 
date for a transition operator. However, positivity and the norm bound ||P|| < 1 
both get lost if we consider the system as a subsystem of a larger one with ob- 
servable algebra, say M n <E> M n . Then (idx„ ®P) takes the unitary flip operator 
$ = Ylij lu)0'*l i n t° n times the one-dimensional projection p = (1/n) 
Hence ||idx„ <g>P\\ > n, and P(l — $) = 1 — np is not positive, although (I — $) is. 
Clearly, id>(„ <S>P is no longer a transition operator, although there is no interaction 
with the "innocent bystander" system described in M. n . 

In order to exclude such phenomena one defines a linear operator P : A — > B 
between C*-algebras to be completely positive [Tak,Pau], if idx n <g>P is positive 
for all n or, equivalently [Tak], if for any choice of n-tuples G A, and 

61, . . . ,b n e B, the operator . b*P(a*aj)bj is positive. P is said to be completely 
bounded, if ||idx„ ®P|| is bounded by a constant independent of n. For such oper- 
ators we define the completely bounded norm as 

||P|| cb :=sup||id M „®P|| . 

If ||P|| cb < 1, P is called a complete contraction. The appearance of the matrix 
algebras M. n in these definitions is solely a matter of convenience: these definitions 
imply the corresponding statements with A4 n replaced by an arbitrary C*-algebra 
M. (see the proof of Lemma 8 for a very similar argument). For checking complete 
positivity or boundedness it is often useful to consider M n <8> A as the *-algebra of 
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n x n-matrices with entries in A. When B is a finite dimensional algebra containing as 
direct summands at most the k x /c-matrices, it suffices to verify complete positivity, 
or to compute ||-P|| cb , in M. n ® A with n = k [Smi]. In particular, all operators 
between finite dimensional C*-algebras are completely bounded. 

Basic examples of completely positive maps are *-homomorphisms, maps of the 
form A i— > V* AV , and all positive maps with either A or B abelian, which includes all 
states. Completely positive operators are completely bounded, and when 1 e A, and 
P is completely positive, we have ||P(1)|| = ||P|| = ||P|| cb [Pau]. The fundamental 
structure theorem for completely positive maps is the Stinespring Dilation Theorem 
[Sti], stating that every completely positive P : A — > B(H) can be decomposed 
in an essentially unique way into P(A) = V*tt(A)V, where tt : A ^ B(fC) is a 
^representation of A on a Hilbert space /C, and V : H — > /C is a bounded operator. 

Basic examples of complete contractions are differences P = P + — P- of com- 
pletely positive maps with \\P + + P_\\ < 1, and multiplication operators A hh> MA 
where ||M|| < 1. Results analogous to the Stinespring dilation are also available 
for completely bounded maps. However, the uniqueness is typically lost. Thus 
any complete contraction P : A — > B{7i) with A 3 1 can be decomposed as 
P(A) = V*tt(A)V2, with tv a ^representation of A, and V\ and V2 isometries. 
Essentially the same statement is that every complete contraction P : A — > B(H) 
can be realized as the off-diagonal corner of a completely positive map [Pau], i.e., 
there is a unit preserving completely positive map P : M2 ®A — > At 2 <8> B{H), such 
that 



Every completely bounded operator P : A — > i3(7i) is a linear combination of 
completely positive ones. If, moreover, P is hermitian (i.e., P{A*) = P(A)*), one 
can find a completely positive P+ with ||P + || cb = ||P|| cb such that P + ± P are 
completely positive [Wit, Pau]. The same statement holds when B(H) is replaced by 
an arbitrary injective C*-algebra [Wit], but fails in general. Since finite dimensional 
algebras are injective, this covers the applications of this result in connection with 
Lemma 14. 
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